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ABSTRACT 


In this research paper, we apply the novel Temimi-Ansari 
method to six first order nonlinear partial differential equations 
for exact solution namely: Burger’s equation, Fisher’s equation, 
Schrodinger equation, wave equation, advection equation and 
KDV_ equations respectively. Unlike other semi-analytical 
iterative methods, this method doesn’t require linearization, 
perturbation, discretization, or the calculation of an Adomian 
polynomials for nonlinear terms in the Adomian decomposition end 
method (ADM). It gives the closed form solution of the problem 
if it exists in finite steps of a converging series that’s 
computationally convenient, easy to obtained and elegant. It 
solves the inherent problem of dealing with the nonlinear term in 
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a straightforward way without stress. The result obtained 


revealed, all the chosen problems give rise to their closed form 
solution in simple steps which confirmed the method is 
powerful, reliable and has wide applicability to other nonlinear 


problems. 
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INTRODUCTION 

Nonlinear partial differential equations technically 
called evolution equations are equations that 
constitute the dissipative term and partial derivative 
of the dependent variable with respect to one or 
more of the independent variables. These equations 
feature prominently in most physical phenomena 
especially in the physical and applied sciences 
(Wazwaz, 2009). They modelled phenomena in the 
field of sciences, engineering, Biology, 
hydrodynamics, chemistry, physics, optical fibre, 
chemical kinetics, plasma physics, medicine and 
may others, Bekir, Tascan and Unsad (2015), 
(Feng, 2002). (Ebiwareme, 2021), Ebiwareme and 
Ndu (2021), (Liberty, 2021). Academics have 
devoted time to extensively studied these equations 
for analytical or approximate solution owing to 
their importance. All though most of the equations 
have closed form solutions, whereas others have 
solution which are difficult to obtained however 
due to the advent of portable computers, the 
solutions can be obtained with less computations. 
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Some of the methods proposed to study these 
equations include: Tanh-Coth method, Asokan and 
Vinodh (2017), (Wazwaz, 2008), (Wazwaz, 2009), 
Chebyshev collocation method Gupta and Saha 
(2015), Vaganan and Asokan (2003), Direct 
similarity method, Manafian, Mahrdad and Bekir 
(2016), G/G' expansion method, Fan and Hongqing 
(1998), Homogenous balance method, (Wazwaz, 
2007),extended tanh method, Fan and Zhao (2000), 
Gomez and Salas (2010), Variational Iteration 
method, Nourazr, Mohsen, Nazari-Golshan (2015), 
Homotopy perturbation method, Hashim, Norami 
and Al-Hadja (2021),Adomian decomposition 
method, Necdet and Konuralp (2006), Darvishi, 
Kheybari and Khani (2008),spectral collocation 
method, Khatter and Temsah (2017), linear 
superposition method, Ma and Fan (2011),trial 
equation method, Gurefe, Sonmezoglu and Misirli 
(2011),Fractional sub-equation method, Zhang and 
Zhang (2011),First Integral Method, Zhang, Zhong, 
Shashan, Liu, Feng and Gao (2013), (Feng, 2013), 
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Ebiwareme and Ndu (2021), Banach contraction 
method, (Ebiwareme, 2021). 


Most recently, Temimi and Ansari proposed a semi- 
analytical iteration method called (TAM) to 
successfully solve to linear and nonlinear functional 
problems. Unlike the ADM, HAM, VIM and HPM, 
this method solved the difficulties that often arise 
for nonlinear terms by writing it in the form of 
Adomian polynomials (ADM), construction of an 
Homotopy (HAM), calculation of the Lagrange 
multiplier (VIM) and tedious calculation of an 
algebraic calculation from corresponding terms as 
in HPM are overcome. TAM has been successfully 
applied to solve various problem, such as ODEs 
Azeez and Weli (2017), Duffing equation, Al- 
Jawary and Al-Razzaq (2016),Chemistry problem, 
Al-Jawary and Raham (2016), Nonlinear ODEs, 
Temimi and Ansari (2015), Thin film flow 
problem, (Al-Jawary, 2017),second order multi- 
point boundary value problems, Al-Jawary, Radhi 
and Ravnik (2017), Fokker-Planck’s equation, 
Temimi and Ansari (2011), Linear and nonlinear 
ODEs, Azeez and Weli (2017),Newell-White-head 
equation, Latif, Salim, Nasreen, Alifah and 
Munirah (2020). 


In this present article, our motivation is to apply 
this novel semi-analytical iteration method to solve 
six differentials nonlinear PDEs. To confirm the 
accuracy, reliability, robustness, and efficiency of 
this method, we seek to ascertain whether a closed 
form or analytical or an approximate solution will 
result from this method. The study is organized as 
follows. In section, the introduction of the study 
detailing PDEs and their applications in physical 
phenomena together with methods used hitherto to 
seek for their solutions. The fundamentals of the 
novel Temimi-Ansari method and the condition for 
its convergence is presented in section 2. In section 
3, we apply the AIM to six different nonlinear 
PDEs and seek closed form solutions and finally the 
conclusion of the study is given in section 4. 


BASICS OF THE ANALYTICAL ITERATION 
METHOD (AIM) 

Consider the general functional 
equation in operator form as follows 


L(u(x)) + N(u(x)) + f(x) =0, (1) 


d , 
B (u =) = 0, or u,(0) = a and u,(0) = b (2) 


differential 


Where x is the independent variable, u(x) is an 
unknown function, f(x) is a given known function, 
L is a linear operator, N is a nonlinear operator and 
B is a boundary operator. 


To implement the TAM method, we first assume 
thatw9(x) is an initial guess that satisfy the problem 
in Eq. (1) subject to Eq. (2). 


L(up(x)) + f(x) = 0,B (uo, 2) = Oor uo(0) = 
aand u (0) = b (3) 


The next approximate solution is obtained by 
solving the problem 
du, 


L(u,(x)) + N(uo(x)) + f(x) =0,B (us) _ 


dx 
0, or u,(0) = aand u,(0) = b (4) 


The next iterate of the problem become 


L(uz(x)) + N(uz(x)) + fx) = 0,8 (w.,“) e 


d 
0, or u2(0) = aand u,(0) = b (5) 


Continuing the same way to obtain the subsequent 
terms, the general equation of the method becomes 


L(tn41(x)) + N(un(x)) na f (x) = 


dunt+1 
0,B (une = 0, or 


Un+1(0) = b 

Then the solution of the problem in Eq. (11) is 
given by 

u(x) = limp. Un (x) (6) 


From Eq. (6), each y(x) is considered alone as a 
solution for Eq. (1). This method easy to 
implement, straightforward and direct. The method 
gives better approximate solution which converges 
to the exact solution with only few members. 


NUMERICAL APPLICATIONS 

In this section, we apply the AIM to solve six 
Partial differential equations that finds usual 
applications in Science and Engineering. They 
include Fishers’ equation, wave equation, advection 


Un4+1(0) =a and 


equation, Korteweg-Devries equation (KDV), 
Burger’s equation and Schrédinger equation 
respectively. 


Example 1. Consider the one-dimensional Burger’s 
equation of the form 


Up + UU, — Uyy = 0 (7) 
Subject to the initial condition 
u(x, 0) = 2x (8) 


Applying TAM to both sides of the equation, we 
get 

L(u) = uz, N(u) = uu, — Uyy, f(x, t) = 0 

The initial problem to be solved is of the form 

L(uo (x, t)) + f(x,t) =0,ug(x%,t) =2x (9) 


Integrating both sides of the above equation subject 
to the initial condition, we get the initial solution as 


Up (x,t) = 2x 
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The second iterative solution is obtained using the 
equation 

L(uy(x, t)) + N(uo (x, t)) + f(x,t) = 

Oxf) = 2K (10) 
Integrating both sides of the above using the initial 
condition yield the integral 


t t 
| via (x, t)dt = | (ton — Ugttog dt 
0 0 


Solving the above integral yield, the second 
iterative solution as 
u,(x,t) = 2x — 4xt (11) 


The next iterate of the problem is given by 

L(u2(x, t)) + N(uz(x, t)) + f(x, t) = 0,u,(%, t) 
= 2x 

Taking the inverse operator of both sides of the 

above yield 


t 
. Uz¢(x, t)dt = J (Uyxx — UyUyx)dt (12) 


Plugging in the derivatives and evaluating, we 
obtain the third iterate as 
6 


Uz(x,t) = 2x — 4xt + 8t?x — 3x (13) 


Similarly, the next iterate is obtained with the 
problem 

L(us(x, t)) + N(u2(x, t)) + f@, t= 

0,u3(x,t) = 2x (14) 


Integrating both sides using the initial conditions, 
we get the fourth iterate as follows 


I U3¢(x, t)dt = FCs — UzU2x)at (15) 


Substituting the derivatives, we obtain the solution 
as 

U3(x,t) = 2x —4xt + Bxt? — 16xt? + =xt* — 
Fe fp ots (16) 

3 9 

Using the relation, u(x,t) = lim,_,.. u,(%, t), the 
closed form solution of the problem is obtained 
u(x,t) = 2x(1 — 2t + 4t? — 8t? +---) 

2x 


u(x,t) = TS (17) 
Example 2. Consider the Fisher’s equation as 
follows 

Up = Ux, + u(1 — u) (18) 


Subject to the initial condition 
u(x,0)=a (19) 
To implement TAM, we have the following 


L(u) = uz, N(u) = —(Uyx aa u(1 _ aL) ef ei t) 
=0 


The first problem to be solved is given by the 
equation 
L(uo (x, t)) + f(x, t) = 0,u9(x%,t) =a (20) 


Taking the inverse operator of both sides yield the 
first iterative solution as 

Ug (x,t) =a 

The second iterative solution is obtained using the 
relation 

L(u,(x, t)) + N(uo (x, t)) + f(x, t) = 

O16 (4t) =e (21) 


Integrating both sides of the above using the initial 
condition yield the integral 


t t 
| Uz(x, t)dt = | (Uoxx — UpUox at 
0 0 


Solving the above integral yield, the second 
iterative solution as 


u,(x,t)=a+a(l-—-a)t (22) 


The third iterate is solved using the relation 

L(up(x, t)) + N(uz(x, t)) + f(x, t) = 0,u,(%, t) 
=a 

Taking the inverse operator of both sides using the 

initial condition, we get 


i Uz(x, t)dt = J, Qs — UjUyx)dt (23) 


Solving the above after plugging in the derivatives, 
we get the solution as 
u,(x,t) =a+tal-—a)t+al—-—a)a(1— 

t? t3 
2a) = + |-a7(1 - a)? o| (24) 
The fourth iterative solution is obtained using the 
problem 
L(us(x, t)) + N(u2(x, t)) + f(x,t) = 
0,u3(x,t) =a (25) 


Integrating both sides using the given initial 
conditions, we obtained the integral as follows 


f, Uz¢(x, t)dt = Gee — UzU2x)at (26) 


Solving the integral above yield the iterative 
solution below. 
u3(x,t)= ata(l-—a)t+a—-—aall— 


2 3 
2a) — + a(1—a)(1—6a+ 6a?) — + ae 07) 


Continuing in similar fashion the subsequent terms 

will be determined, hence the solution of the 

problem is obtained using the relation 

u(x,t) = lim u,(%,t) 

n—o0 

t2 

u(x,t) = a+a(1—a)t+a(1—a)a(1 - 2a) Tt 
t3 

+a(1—a)(-—6a+ 6a*) + ee 
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aet 


u(x,t) = 


(28) 


1-ataet 


Example 3. _Let’s consider the first-order wave 
equation in one-dimension as follows 
uz + ku,, =0,k > 0 (29) 


Subject to the initial condition 

(x, 0) = sin (=) u,(0,t) = are (=) 
BAP adele sere a 

: —krt -k x 

u(0,t) = sin (—) ue, 0) = = COs (=) 

(30) 


Applying TAM to both sides, we get the following 
terms 


TCX 
L(u) = ur, N(u) = Kuyx, Up (x, 0) = sin (—) 


The first iterate is obtained by solving the problem 
of the form 


L (u(x, t)) + f Ge, £) = 0, uo x,t) = sin (™) 
(31) 


Integrating both sides subject to the initial 


condition, we get the initial solution as 


U(x, t) = sin (=) 


The second iterative solution is obtained by solving 
the problem 

L(u(x, t)) + N(uo (x, t)) + f(x, t) = 

0, u, (x,t) = sin (=) (32) 


Integrating both sides of the above using the initial 
condition yield the integral below 


I Uz,(x, t)dt = f Gos — UpUox)at (33) 


Solving the above integral yield, the second 
iterative solution as 


wea = [sn()-cos(@))(#) 

The next iterate of the problem is given by 

L(u2(x, t)) + N(uz(x, )) + f(x, t) = 0,u2(x, t) 
x 


=sin(% 


Taking the inverse operator of both sides of the 
above yield 


i Uz(x, t)dt = cre — Uz Uy x)at (35) 


Plugging in the derivatives and evaluating, we 

obtain the third iterate as 

Uz(x,t) = 

. TX 1 (knt\? mx \ | krt 1 (knt\3 

sn (2) [1-3 2)’ ]-om 2) EES 

l aX l Path 

(36) 

Continuing in the same way, the succeeding terms 


are obtained, and the closed form solution of the 
problem is given by 


u(x,t) = sin (=) cos (=) 


Using trigonometric identity, the above reduced to 
the form 


u(x,t) = sin [x (—)| (37) 
Example 4. Let’s consider the homogenous 
advection equation of the form 

Ut + UU, = 0 (38) 


Subject to the initial condition 

u(x,0) = —x (39) 
Implementing TAM on both sides of the equation, 
we have the following terms 

L(u) = uz, N(u) = uu,, f (x,t) = 0 

The first problem to be solved for the first iterate is 
given by 

L(uo (x, t)) + f(x,t) = 0,ug(x,t) =—x (40) 


Integrating both sides of the above equation subject 
to the initial condition, we get the initial solution as 
Ug (x, tp=Sx (41) 
The second iterative solution is obtained using the 
equation 

L(u(x, t)) + N(uo (x, t)) + f(x,t) = 

0,tij (x, t) = -—x (42) 
Integrating both sides of the above using the initial 
condition yield the integral equation 


t t 
| U44(x, thdt = -| (uUgUpx)at 
0 0 


Solving the above integral yield, the second 
iterative solution as 
u,(x,t) = —x — xt (43) 


The next iterate of the problem is given by 


L(u2(x, t)) + N(uz(x, t)) + f(x, t) = 0,u,(%, t) 
=—x 


Taking the inverse operator of both sides of the 
above yield 

t t 
J Ug¢(x, t)dt = — J (UU; ,)dt (44) 
Plugging in the derivatives and evaluating, we 
obtain the third iterate as 
(4,0) S=—X = xt = x07 = =xt3 (45) 


Similarly, the next iterate is obtained by solving the 
problem below 


L(ug(x, t)) + N(ug(x, t)) + f(x, t) = 0, u3(x, t) 
= —x 


Integrating both sides using the initial conditions, 
we get the fourth iterate as follows 
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t t 

Jy Use(x, dt = — J, (uguzx)dt (46) 

Substituting the derivatives, we obtain the solution 

as 

U3(x,t) == —x — xt — xt? — =xt3 — =xt4 —e 
(47) 

Using the relation, u(x,t) = lim,_,.. u,(%, t), the 

closed form solution of the problem is obtained 

u(x,t) =—x(1+t+t?7 +t? 4+t*+--) 

u(x,t) = (48) 


t-1 
Example 5. Consider the Korteweg-de Vries 
(KDV) equation which takes the form 
Ut — GU, + Uyxy, = 0 (49) 


Subject to the initial condition 

_ ke 2 (kx 
u(x, 0) = —7sech (=) (50) 
Applying TAM to both sides of the equation yield 


the expressions 
L(u) = uz, N(u) = —6uu, + Uzyy, f (x, 2-0 


The initial problem to be solved is of the form 
L(uo(x, t)) + f(x, t) = 0,u9(x,t) = 

k? 2 (kx 
—7sech (=) (51) 
Integrating both sides of the above equation subject 
to the initial condition, we get the initial solution as 


U(x, t) = —““sec h2 (=) (52) 


The second iterative solution is obtained using the 
relation 


L(u,(x, t)) + N(uo (x, t)) + f(x, t) = 
ee 2 [kx 
0, u(x,t) = —sech (=) (53) 
Integrating both sides of the above using the initial 
condition yield the integral 
t t 
J U(x, t)dt = J (6UgUox — Uoxxx)dt — (54) 


Solving the above integral yield, the second 
iterative solution as 


u,(x,t) = 
5 2(kx 
— sec h? (2) — at tan h (2) (55) 


The next iterate of the problem is given by 


L(u2(x, t)) + N(u(x,t)) + ft) = 
Otol x,t) = 2% (56) 


Taking the inverse operator of both sides of the 
above yield 

t t 
J Ug, (x, t)dt = So (6u, Uy — Urxxx)at (57) 
Plugging in the derivatives and evaluating, we 
obtain the third iterate as 


U2 ‘ea t) = 

k2 2 (kx k5 sech?(“*) eae 
see (t) —# PE ay (8) 
—sech* |= (2 — cosh[kx])t? + = (58) 


Continuing in the same, the succeeding term are 
obtained in a similar manner 


Using the relation, u(x,t) = lim,_,.. u,(%, t), the 
closed form solution of the problem is obtained 


tint) = —"sech? [F (x — k?t)| (59) 


Example 6. 
in the form 
uz + tux, = 0,u(x,0) = 14+ cos h(2x) (60) 


Consider the Schrédinger equation 


Applying TAM to both sides of the equation, we 
get the constants 

L(u) = uz, N(u) = iuy,, f(x, t) = 0 

The first problem to be solved is given by 

L(up (x, t)) + f(x,t) = 0,ug(x,t) =1+ 

cos h(2x) (61) 


Taking the inverse operator of both sides subject to 

the initial condition, we get the initial solution as 

Uo (x,t) = 1+ cos h(2x) (62) 

The next iterative solution is obtained using the 

equation 

L(uy(x, t)) + N(uo (x, t)) + f(x, t) = 0,u,(x, t) 
= 1+cosh(2x) 


Integrating both sides of the above using the initial 
condition yield the integral 


t Ce 
Jy Uae(x, thdt = J, (ittoxx)dt (63) 
Solving the above integral yield, the second 


iterative solution as 
u(x,t) =1+cosh(2x)+4icosh(2x) (64) 
The third iterative solution of the problem is given 
by 
L(u2(x, t)) + N(uz(x, t)) + f(x, t) = 
0, u(x,t) = 1+ cos h(2x) (65) 
Taking the inverse operator of both sides of the 
above yield 

t f 
| Up, (x, t)dt =| Cite Jat 

0 0 


Plugging in the derivatives and evaluating, we 
obtain the third iterate as 
Uuz(x,t) = 1+ cosh(2x) + 4icosh (2x) + 
aey2 
(4it)" cos h(2x) (66) 


2! 


Similarly, the next iterate is obtained with the 
problem 
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L(ug3(x, t)) + N(ug(x,t)) +f, 0) = 
0,u3(x,t) = 1+ cos h(2x) (67) 


Integrating both sides using the initial conditions, 
we get the fourth iterate as follows 


t Ce 
Sy Uz¢(x, t)dt = Sy (iugxx)adt 
Substituting the derivatives, we obtain the solution 


as 
u3(x,t) = 1+ cosh(2x) + 4icosh (2x) + 


. 2 . 3 
— cos h(2x) + — cos h(2x) 


(68) 


(69) 


Using the relation, u(x,t) = lim,_,.. uy (x, t), the 
closed form solution of the problem is obtainedas 


_» , (it)? 
u(x,t) = 1+ cosh(2x){ 1+ (4it) + 1 
(4it)? (4it)* ) 
u(x, t) = [1 + cos h(2x)]e** (70) 
CONCLUDING REMARKS 


In this research article, the closed form, or exact 
solutions for six different nonlinear partial 
differential equation is investigated using the novel 
Analytical Iteration Method (AIM). The efficiency 
of the method is confirmed by solving the Burgers, 
Fisher’s, Advection, Schrddinger, wave and KDV 
equations respectively. The method gives an 
analytical solution which converges rapidly to the 
exact solution with more terms considered. This 
solution is easily verifiable in few steps of iteration 
subject to the initial condition and _ is 
computationally convenient since it doesn’t require 
perturbation, discretization, and linearization. It is 
observed the method is reliable, efficient, and 
applicable to all class of nonlinear problems in the 
fields of Science and Engineering. 
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